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The ratio of the a-decay widths of a bare nucleus and the related atom is calculated. Both the 
change of the form and thus the penetrability of the potential barrier and the effect of reflection 
in the classically-allowed region appearing due to the electron shell are taken into account in the 
calculations of this ratio. The contribution of each of these two effects is of one and the same order of 
magnitude. For long-lived radioactive samples the values of the total effect turn out to be somewhat 
below 1 percent. 

PACS numbers: 23.60.+e 



I. INTRODUCTION 

Transmutation of long-lived radioactive isotopes is one 
of the hottest problem of the modern nuclear physics. Re- 
cent attempts to produce an affect upon the decay rate 
of such isotopes by physical methods which do not deal 
with any beam of neutrons, gamma-quanta or charged 
particles demonstrate increasing interest and popularity 
of this problem. By now a list (not so wide) of exam- 
ples demonstrating the possibilities to vary the rate of 
/3-transitions by deep ionization of atoms is known 0, Q ■ 
Also it is shown that the half-life times of some low- 
energy 7-transitions turn out to be much smaller when 
the energy of the nuclear transition is close to the energy 
of any allowed atomic transition of the same multipolar- 
ity [3j, |j| . The a-decay looks as a process which is the 
most insensitive to the external influence of the discussed 
type due to the typical value of the energy of the potential 
barrier which is in the neighborhood of 20 MeV. That is 
why this process attracts the smaller interest than other 
types of decay. Nevertheless in the last decade some the- 
oretical papers devoted to this problem were published 
[H, Moreover an interesting supporting evidence is 
obtained in the experiment presented in Ref. [7]. The 
difference between the decay widths of 212 Po in Pb and 
Ni matrixes is proved there with 95 per cent reliability 
(see the discussion bellow). Thus, in our opinion, it is 
necessary to study accurately the effects modifying the 
a-decay widths. 

In the present paper the a-decay widths of a bare nu- 
cleus and the respective atom are compared. The elec- 
tron cloud changes, first, the kinetic energy of the emitted 
a-particle and, second, the decay rate. The effect of the 
kinetic energy difference in cases of the decays of the bare 
nucleus and the atom is well-known from the middle of 
the last century. As usual nuclear data tables contain the 
values of energy yield of the a-decay of atoms. Therefore 
the effect can not be neglected and should be accurately 
considered: the formal inclusion of such energy correc- 
tion leads to increasing of the a-decay width by the fac- 
tor lying between 1.22 (for 212 Po) and 2.60 (for 147 Sm) 
Q. This circumstance is always taken into account. The 
problem how the decay rate is changed by the electron 



cloud has not solved accurately up to now. 

A closely related problem is examined in Ref. [6]. To 
some extent paper [6j is a reaction of scientific commu- 
nity to mass media publications in which the a-decay 
enhancement is anticipated in the conditions, similar to 
that discussed by the proponents of the so-called "cold fu- 
sion". The paper provides a qualitative argument against 
these anticipations. It is shown that the upper limit of 
the ratio of the a-decay half-lives of a nucleus screened 
by electrons of a surrounding dense medium and the re- 
spective unscreened nucleus in actual cases is not larger 
than 1.02. At the same time the author does not de- 
cide between two results obtained for two versions of the 
screening (~ 1.013 and 1.000 in 226 Ra) i. e. between the 
modest and zero effects. Some other details of the in- 
vestigations which is carried out in Ref. |(| are discussed 
bellow. So in our opinion the outlined above problem has 
been unresolved till now. 



II. STATEMENT OF THE PROBLEM AND 
FORMALISM 

It is convenient to begin consideration of the prob- 
lem with a simple quasi-classical approximation. The 
potential of the interaction of the a-particle with the 
atomic shells is monotonically increasing negative func- 
tion V s h(r), very slowly varying in the potential barrier 
region. Neglecting the r-dependence of the potential in 
this region (assuming V s h(r) = — Vo), one obtains the 
only effect of the shift (equal to Vo) of the decay energy. 
As it is pointed out in Ref. Q the a-decay width re- 
mains the same. A small increase of the V s h (r) results in 
increasing of the barrier height and consequently in the 
decreasing of the barrier penetrability which takes the 
form: 



P at = exp 



Text 

-~ J VVbn (r) + V sh (r)- E bn + V dr 



, (1) 



where Eb n is the resonance energy in the bare nucleus 
case, Vb n [r) = V str (r) + V cf + V Co ui{r) is the sum of the 
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strong, the centrifugal, and the Coulomb potentials of the 
interaction of the a-particle and the residual nucleus. 

At the same time there is another effect originated 
by the atomic shell. It appears due to the fact that 
the asymptotic behavior of the resonance (Gamow) wave 
function in the case of bare nucleus differs from that of 
the respective atom. Indeed, the asymptotic behavior of 
a resonant function looks as follows: 

X[ es (r) =D[Gi {i],kr)+iFi {r),kr)], r -> oo, (2) 

where Gi(r/,p) and Fi{rj,p) are the Coulomb wave func- 
tions (Tlj . 



V = (e 2 /hc)Z 1 Z 2V / pc 2 /(2E) 



(3) 



is the Coulomb parameter, eo > - the elementary 
charge, E - the decay energy yield, k - the respective 
wave number, and p - the reduced mass, Z 2 = 2 - the 
a-particle charge. The coefficient 



(4) 



D = y/flT/h 



determines the decay width T. In the case of bare nu- 
cleus the parameters Df, n , rjbn, and kf, n are related to 
the energy Eb n and the charge of the daughter nucleus 
Z\ = Z—2 whereas the corresponding parameters rj at and 
k at (the case of neutral atom) are related to the energy 
E a t = Ei, n — Vo and the charge Z\ = — 2 of the residual 
system - the nucleus and the electron cloud. The origin 
of the discussed effect is the reflection of the a-particle 
wave on the varying potential. An obvious way to ac- 
count for this effect is to solve the Schrodinger equation 
directly. It may be found grater or smaller than the ef- 
fect of the barrier form variation considered above. The 
main goal of the present paper is to calculate the ratio of 
the decay width, accounting for both these effects. 

A physical model applied here for the description of 
the a-decay of the atomic system includes the following 
ingredients. 

First, the distribution of the electronic cloud is as- 
sumed to be undistorted during the a-particle emission. 
This assumption is based on the facts that the orbits of 
the strongly bound electrons are slightly changed due to 
the variation of the nuclear charge Z —> Z — 2 and the 
weakly bound electrons are too slow. Indeed, the elec- 
tron velocity v e rs v a (v a is a velocity typical for the a- 
particle) corresponds to the electron energy E e ~ 500 eV. 
Thus a high speed and a small charge of the emitted par- 
ticle enables one to represent the potential energy term 
of the Hamiltonian in the folding-form interaction of the 
a-particle with negative ion consisting of the daughter 
nucleus and the atomic shell of the mother one. 

Second, the electron distribution is taken in the 
Thomas- Fermi form [9|. 

Third, the electron charge, associated with sphere of 
the radius Rq which is somewhat greater than that of the 
range of strong interaction R str ~ [1.2 J 4 1 / 3 +i? Q + 2] fm, 
is assumed to be small enough to neglect the variation of 
the potential originated by this charge. 



At last only the ratio of the a-decay widths but not the 
absolute values of the widths for both an atom and a bare 
nucleus is the object of the calculation. Evidently the 
absolute values are related to many-nucleon properties of 
a nucleus and thus are beyond the scope of the present 
paper. 

The properties of the Gamov solutions to the two-body 
Schrodinger equation related to the presented model 
are the following. The solution is characterized by the 
asymptotic form @ with the parameters D atl rj a t, and 
k at at very large (several times grater than the atomic ra- 
dius) distances. At smaller distances it can be obtained 
only numerically. In the region R str < r < Ro the equa- 
tion turns out to be the Coulomb one with the parameters 
rjbn and kbn related to the bare nucleus. Consequently in 
this domain the resonance solution becomes: 

X[ es (r) = Dot [AG i {rjbn , k bn r) + BFi (r? 6 „ , k bn r)] . (5) 

In this region the coefficients D a t, A, and B differ the 
solution from the one, corresponds to the Gamow state 
of bare nucleus, where D = Db n , A = 1, and B — i (see 
Eq. ©). 

In the both cases the amplitude of each wave func- 
tion is determined by the matching of the function of the 
type ([5]) with one and the same internal wave function. 
As it is shown in the Ref. the term containing the 
regular Coulomb function plays a negligible role in the 
decay width calculations of a narrow resonance. Indeed, 
the ratio of the regular and the irregular Coulomb wave 
functions in the sub-barrier region can be estimated as 
follows: 

FiiVbn, h n r) /Gi(T] bn , k bn r) = P(r) = 



exp 



y/E - V{p) dp 



r < Text- 



(G) 



Remind that for the a-emitters with the half-lives r > 
1 y the value of the total penetrability of the barrier 
(which is expressed by Eq. © with the lower limit r = 
Tint) P < 10~ 27 . Furthermore, the region beyond the 
radius R s t r contributes predominantly to the value of the 
integral in Eq. ©• Thus, the second term in Eq. ([5]) 
can be neglected with a high precision. Therefore the 
squared ratio of the amplitudes of the wave functions 
and consequently the ratio of the widths become 



r 6 „/r Qt = \D bn /D at \ 2 = \a\ 2 



(7) 



This coefficient should be calculated numerically. The 
calculation procedure is complicated due to the need to 
solve numerically Schrodinger equation in extremely wide 
domain of values of the radial variable r. The results of 
the performed calculations demonstrate that the asymp- 
totic form © becomes valid at the distance of the order 
of 10 6 fm whereas a typical value of the wave number of 
the a-particle is about k ~ 1 fm . At the same time a 
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high numerical precision should be provided both at the 
large distance and in the neighborhood of the nucleus. In 
fact one deals with two very different scales of distances 
in one and the same calculation procedure. 

The optimal way is as follows. As it is noted above 
one has no need to solve the Schrodinger equation in the 
r < R s tr region for the calculation of the ratio of the 
decay widths of an atom and the respective bare nucleus 
- this ratio is determined beyond the region. Therefore 
it is convenient to consider the reduced equation: 



Gi( Vat ,k at r) = [-7x[ F) (r) +a^( r )]/A 



IG), 



+ Vcoui(r) + V sh {r) + V cf (r) -E) X i(r)=0 



2[i dr 2 

(8) 

The Thomas-Fermi model was utilized for description 
of the electron cloud potential: 



V sh (r) = 2tp sh {r)e , 



Vsh{r) 



(!-£(*)), 



(9) 



(10) 



x = rZ 1/2 m e elb~ 1 h~ 2 : , b = (3tt/4) 2/3 /2 (m e - electron 
mass). In this case Vo = 2m e Z 4 ' 3 e^Kb^hr 2 and k « 
1.59 is the constant in the asymptotic formula £(x) ~ 
1 — kx ) as x — > 0. We introduce the distance Rq in such a 
way that in the domain r < Rq the r-dependence of V s h 
can be neglected: V s h(r) = —V . Thus in the domain 
R s tr < r < Rq equation ([8]) takes the form 



d 2 xi(r) , {2(,E bn ^(Z-2)e 2 1(1 + 1) 



dr 2 



V h 2 



rr- 



Xi(r) = 0. 

This equation coincides with that in the case of bare nu- 
cleus but with shifted energy Ef, n . 

The Coulomb wave functions Fi(rjbn, kbnr) and 
Gi(r]bn, kb n r) are the solutions to Eq. (|TTj) . Thus the two 
linearly independent solutions to Eq. ([5]) can be chosen 
in the following way 



x[ F '(r) = Fi(rj bn , k bn r), R s tr < T < Rq 



X\ (r) = Gi{ribn, k bn r), R s tr < r < R . 



(12) 



(13) 



At large distances, far from the electron shells the sum 
of potentials of electron shells and nucleus looks like the 
potential of negative ion with the charge Z\ = — 2. In 
this region the solutions to Eq. (jSJ) have the form of linear 
combinations of functions Fi(rj a t, ka.tr) and Gi(rj a t, k a tr): 

xP(r) = aF^ukatr) + pG^aukatr), (14) 



x\ '(r) =jFi(rj at ,k at r) + 5Gi(-q at ,k at r). 
The reversed relationships 

F^uukatr) = [8xf\r)-Pxf\r)]l^ 



(15) 



(16) 



A = aS - (3 j 



(17) 



(18) 



follow directly from Eqs. {H} and (|T5[). The value of A 
can be found without explicit values of a. j3. S, 7. Indeed, 
the relationship 

Wr{xl F \xi G) ) = WriFtiriaukatr^Giiriaukatr^A 

(19) 

follows directly from ([ri]). ITS)) , and fig]). The lower in- 
dex r means that the derivatives with respect to variable 
r appear in the Wronskian W r . As it is known [ll[ 



W p (F l (r,,p),G l (r,,p))=l ) 
so from (fT2]l. ([T3"]). and (O it follows that 

A = hn/kat- 



(20) 



(21) 



In the case of a-decay of a neutral atom the resonant wave 
function has the asymptotic form D a t [Gi(r] a t, k a tf) + 
iFi(r?at, k at r)] at large distances far from the electron 
shells. Now, with the aid of (|16|) . (jTTJ) it is clear that 
the considered resonance wave function looks as follows 



Xi 



D n 



a 



i/3 (G) 
— Xi 



.S + vy (f) 
1 ^^Xi 



(22) 



The function (|22p in the domain R str < r < Rq has the 
form (see (SJ, (SJ): 



D„ 



— ^ — Gi(T]bn, k hn r) + 1 . Fi(T] bn , k hn r) 



(23) 



Thus, we come to the following conclusion. In the 
case of the neutral atom the internal wave function must 
be matched with the function ([23)) . while in the case of 
the bare nucleus with the function D a t \Gi(ribn, kbnf) + 
\Fi(rjbm kbnf)\ ■ The matching procedure can be realized 
at any point of the domain R s t r < r < R$. 

In the case of neutral atom the multiplier of 
Gi(rjb n ,kb n T) in Eq. (|23p . which is symbolized as A in 
Eqs. ([5]) and (JT]), changes the width of the resonant state. 
Obviously there is no need to solve the resonant problem 
as a whole to determine the coefficient A. It is sufficient 
to solve (numerically) the Schrodinger equation ((HJ in the 
domain r > R s t r with the boundary condition (|12[) . Then 
the coefficients a and f3 can be deduced from the Wron- 
skians of the solution and the functions Gi(r) at , k at r) and 
Fi(j] at ,k a tr). The explicit expression 



r&n/r Q 



la- 



i/3)/A| 2 = (a 2 +/3 2 )/A 2 



(24) 



for the suppression factor - the ratio of the decay width 
of the bare nucleus and the neutral atom - follows di- 
rectly from and §5§ . It should be pointed out that we 
arrived at a surprising conclusion. Indeed, one needs to 
find the solution with boundary condition (|12[) (regular 
function Fi(rj at ,k a tr)) to calculate the coefficient in the 
term containing the irregular function Gi(rj a t,k a tr). 
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III. RESULTS AND DISCUSSION 

The influence of the electron surrounding onto the a- 
decay width was investigated for some isotopes in the 
framework of the rigorous approach described above. 
The Runge-Kutta and Stoermer methods were applied 

Table I: The values of the coefficient of suppression of the 
a-decay widths by the electron shell. 



Isotope 


Qat 




Pbn/ Pat 


226 Ra 


4.7806 


1.0021 


1.006 


232 rj 1 ^ 


4.0288 


1.0052 


1.010 


144 Nd 


1.9251 


1.0079 


1.014 


148 Sm 


1.9858 


1.0085 


1.015 



to solve numerically the radial Schrodinger equation ([H]) . 
The results obtained by means of these methods coincide 
with high accuracy. The factors of the a-decay suppres- 
sion caused by the electron shells are presented in Table 
I. The third column of the table demonstrates the re- 
sults of precise calculations and the last one the results 
of quasiclassical calculations of the barrier penetrabili- 
ties. The results of exact calculations, based on Eq. 
differ from that of quasiclassical approach. The reason 
of this discrepancy is the following. The ratio Pbn/ Pat 
is the exponential function of difference of two integrals 
(see ([T])) over the subbarrier regions. It turns out that 
the vicinity of the corresponding external turning points 
(upper limits of the integrals) contributes predominantly 
to the difference of the quasiclassical integrals. But the 
vicinity of a turning point is just the region, where the 
quasiclassical approximation looses accuracy. 

In a qualitative sense the results presented in the ta- 
ble are not in contradiction with the experimental one 



where the following difference between the half-life 
times of 212 Po decay in Pb and Ni matrixes is obtained: 
T 1/2 (Pb) -T 1/2 (Ni) = (-0.66 ±0.25) ns. This difference, 
related to the half-life time of 212 Po isotope Tj/2 = 0.3 
fj,s is about 0.2 per cent. It should be stressed that the 
settings of the problems in the present work and in Ref. 
Q are different. The effect of the change of solid state 
surrounding of the emitter but not the effect of the ioniza- 
tion is studied in the experiment. So, a direct comparison 
is impossible. 

As it is noted above the upper limit of an a-decay half- 
life time variation in a dense (solid or liquid) medium is 
estimated theoretically in Ref. [H| ■ The effect is assumed 
to be originated by the change of the quasiclassical pene- 
trability of the barrier by the potential of the interaction 
of the a-particle with electrons of the medium. The ini- 
tial barrier is chosen in the simple form of "cut Coulomb 
potential". The change of the barrier is described by 
the exponential screening factor exp(—r/R). Two pa- 
rameters R = 3.7 • 10 3 fm (for so-called "Thomas-Fermi 
screening") and R = 4 • 10 5 fm (for "Debye screening") 
are considered. The values Pbn/Pat = 1-013 and 1.000 
for 226 Ra are obtained in these cases respectively. The 
larger value of the upper limit is qualitatively confirmed 
in our calculations. 

So the ratio of the a-decay widths of a bare nucleus and 
the related atom is calculated in accurate way. The math- 
ematics providing a possibility to take into account the 
effect of reflection of the a-particle wave at the electron 
shell in the classically-allowed region is built. The change 
of the penetrability of the potential barrier is taken into 
account simultaneously For long-lived radioactive sam- 
ples the values of the effect turn out to be somewhat 
below 1 percent. They decrease with increasing of the 
decay energy. 

The authors are thankful to V. G. Kalinnikov for mo- 
tivating discussion. 
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